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Abstract 

We show that parity is conserved in vector-like supersymmetric 
theories, such as supersymmetric QCD with massive quarks with no cu- 
bic couplings among chiral multiplets, based on fermionic path-integrals, 
originally developed by Vafa and Witten. We also look into the effect of 
supersymmetric breaking through gluino masses, and see that the parity- 
conservation is intact also in this case. Our conclusion is valid, when only 
bosonic parity-breaking observable terms are considered in path-integrals 
like the original Vafa- Witten formulation. 
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1. Introduction 



Non-perturbative chiral symmetry breaking [1] is an important aspect for phenomenolog- 
ical model building based on vector-like theories, such as in composite particle models [2] [3]. 
It has been proven that parity symmetry is conserved in (non-supersymmetric) vector-like 
QCD theories, and is not broken spontaneously even non-perturbatively [4] [5]. This proof 
[4] is based on the evaluation of fermionic path- integral, which gives always non-negative 
vacuum energy after adding parity-breaking terms. However, the question has risen whether 
the Vafa-Witten constraint [4] [5] for the non-supersymmetric case can be avoided in super- 
symmetric vector-like theories because of the new interactions among gluino-quark-squarks 
[6] [7] , and whether parity is broken like other global symmetries [8] . These particular inter- 
action terms with scalar-dependence seem to be the main obstruction for the proof for the 
positive defmiteness of the determinant in the fermionic path-integral [5] [6] [7]. Moreover, 
the results about gauge symmetry breaking for massless supersymmetric QCD, when the 
number of flavor Nf is smaller than the number of colors: Nf < N [9] also seem to sug- 
gest the parity-breaking in supersymmetric vector-like theories. Although these observations 
seem reasonable, it also seems to contradict with the other universal wisdom about super- 
symmetry that supersymmetric vacuum is stable when the Witten index Tr (— 1) F [10] 
is non-zero, e.g., Tr (— 1) F = N for the gauge group SU(N), and therefore the vacuum 
energy most probably stays zero with no parity breaking. It has been also recently point out 
[11] that supersymmetric QCD has condenstate-free phase with no gluino condenstate. In 
our present paper we reconsider this subtle problem of parity breaking, and give a proof for 
parity conservation in supersymmetric vector-like theories with massive quarks. 

The model we deal with in this paper is a globally supersymmetric vector-like theory 
with massive chiral multiplets coupled to a non-Abelian vector multiplet with no cubic 
coupling among chiral multiplets. Our proof is based on three major assumptions: First 
one for the massiveness of all the quark chiral multiplets, the second one about the absence 
of the Yukawa-couplings among chiral multiplets, and the third one that we rely on the 
method in [4] for purely bosonic parity-breaking observables. Therefore our method does 
not cover the fermionic parity-breaking observables like Wilson fermions [12] treted in lattice 
QCD.Q The massiveness of quarks are also important for non-perturbative conservation of 
supersymmetry, due to well-defined non-zero Witten index Tr(— 1) F in such cases [10]. 
This is because supersymmetry is conserved, only if the vacuum energy is zero. Therefore the 
non-perturbative breaking of supersymmetry would cause the shift of vacuum energy, causing 
the breaking of parity symmetry [4] [5]. Interestingly, we will find that parity is conserved 
also for a vector-like supersymmetric theory, like non-supersymmetric vector-like theory. We 
also look into the effect of gluino masses, which will not disturb the main body of the proof 
for supersymmetric case, and therefore parity is also conserved in broken supersymmetric 
vector- like theories. 



2 For reviews for supersymmetric QCD on lattice, see, e.g., [13] 

2 



2. Review for Non-Supersymmetric Vector-Like Theory 

We start with reviewing the parity conservation in non-supersymmetric case [4] [5] first, 
in order also to elucidate our notation. Suppose the total lagrangian C(X) = C — XX with a 
parameter A is a generalization of the parity-conserving lagrangian C, such as that of QCD, 
with a parity-non-conserving observable X, such as the FF -term, with a real constant A. 
If parity is broken in the vacuum and (X) ^ 0, then the theory can choose a vacuum state 
in which X(X) < due to the signature ambiguity of (X). Hence the vacuum energy 
E(X) can be lower than E(0) of the parity-conserving vacuum energy: E(X ^ 0) < E(0). 

However, an explicit evaluation of path-integral reveals that this would not happen, i.e., 
there is no such vacuum whose energy is lower than that of the parity-conserving one [4]. 
Consider the path-integral in Euclidian space for the vacuum energy E(X): 



e -VE(x) = J [d^/] [ dx ] [dx] [df] [d£] exp [- J d'x (C + iXX) 



(2-1) 



where V is the Euclidian volume, Aj is the gluon field. The indices /, j, - are 
for adjoint representations of the gauge group G. For example for G = SU(N), we have 
/, j, ■■■ = i, 2, n 2 -i. In order to clarify basic constituents of our system, we use 2-component 
spinors in this paper. Since we are dealing with a vector-like theory, our two-component Weyl 
spinors X- an d € with the flavor indices i, j, ■■■ = l, 2, n for quarks are in the conjugate 
representations to each other. The factor of i in the AX -term is due to the usual Wick 
rotation. We specify the lagrangian as 

£ = — \iFab 1 ) 2 + £>F , 

C F = + i(^ i p^x pi )+Kr i pJ^^+m^(x^^ j )+mi(C j xJ , (2.2) 

where m = (m, J ) is an N x N hermitian mass matrix: (r?V)* = to/, which can be 
arranged to have only positive eigenvalues. Since we are dealing in this paper only with a 
vector- like theory, the quark fermions X- an d C are in the representations conjugate to each 
other, e.g., N and N* -representations of SU(N) , respectively. Accordingly, our covariant 
derivative D a contains the minimal coupling of the gauge field to these fermions. We are 
using the notations similar to that in [14], e.g., we use the Minkowskian four- dimensional 
(4D) vector indices a, b, ■■■ = 0, 1, 2, 3, with the signature (i] ab ) = diag. (+,—,—,—), while 
a, p, ■■■ = 1, 2 and 0, ■■■ = 1, 2 for the 2-component spinors. Other relevant relations are 
such as 



'a/3 



(3a ' V VP ' ^a ^ ^p' 



dl = ~d a , (Cap) = (C..) = ( ° 7^ . (2.3) 



Based on these, it is easy to confirm the reality of each term in our lagrangian. 

As was shown by Vafa-Witten in a vector-like theory [5], the fermionic space for the 
path- integral (2.1) can be a direct sum of the positive and negative eigenstates of the Dirac 
operator in a finite volume V. Let x ai an d C^ ifi correspond to such eigenstates [15]: 

We then easily see that the set of X a ^_^ = X^+^y iai{-ii) ^ corresponds to the 

eigenvalue — /x: 

Therefore the whole fermionic space is not only a direct sum of positive and negative eigen- 
states, but also they are always paired up between + /i > and — /x < 0. 

Or equivalently, in terms of a four-component Dirac spinor 

^ = ^ai)=(^) , ^=(^)=(^] , (2.6) 
with the four-component spinorial indices a = (a, a), p = (p,p), -, we have 

V)© ■ MA V)- 

We can also introduce the 2N x 2N mass matrix m and the usual 4x4 75 -matrix for 
the four-component notation by 

~ fm i i \ (5J* \ . . 

Note also that m is hermitian: = m. The eigenstates in (2.6) are much transparent 
now for p as 

Pipy. = +M» > ^V'(-Ai) = -/^(-m) > ( 2 - 9 ) 

The eigenstates V'(-m) correspond to V'(-m) = 7sVyj because 75 satisfies {75,^} = 0, and 
therefore 

IV(-ji) = P(l5^) = ~li>P% = -^75^ = • (2-10) 

Therefore the pairing between /x and — /x eigenstates in (2.6) is clear. Accordingly, the 
lagrangian C-p is simply 

£ F = ^(i|) +m)ij , (2.11) 
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Our path-integral (2.1) is now 

e- yE ^ =J[dA a I ] J[d^][d^]e'I d4x ^ +c ^ xx ) 
= J [dAj] J [#][#] e -Jd*x(c B+ i\x) exp 

= J [dAj] I F exp - J d 4 x(C B + iXX) . 
Here I F is the fermionic determinant from the fermionic path-integral 



d 4 xip(ip + m)ip 



(2.11) 



I F = J [dip] [dip] exp 

= Det (ip +m) 

= Y[ det (ipLhN + m) 

= (detm) nL0+nR0 



J d A xip(ip +m)ip 



[J det (i(J,I 2 N + m ) 

At>0 



\\ det (ifiI 2 N + m) 



At<0 



\2n 



\~[ det (ifxhN + m) 

At>0 



= (detm 

= (detm) 2n ° TT det (ifiI 2 N + ^) det (— i^xI 2 N + ^) 



Q det (- i///2jv + rn) 

M>0 



At>0 



= (detm) 2n ° Y[ det (ifJ,I 2 N + m ) det (+ ifJ,I 2 N + m)* 
= (detm) 2n ° ]~] [det (i/iI 2 N + m) } [det (+ i/iI 2 N + ra) 



(2.12a) 



(2.12b) 



= (detm) 



2n 



J"J det (i/jLl 2 N + m) 

M>0 



> 



(2.12c) 



Here I 2 n is an 2N x 2iV unit matrix, and in (2.12a), no is the number of u = modes, 
satisfying n 



LO 



n R0 = uq. This is because we have to consider only the instanton number 



zero background n L(j — n RQ = that is connected with the original vacuum with energy 
E(X = 0) [16]. The previously-mentioned parings u <-> — u are used also in (2.12a). The 
determinants in (2.12a) are taken for the 2N x 2N matrix for flavour indices, distinguished 
from the symbol 'Det' for the fermionic path-integral. In (2.12b) we have also used the 
hermiticity of m. Now the positive definiteness of I F is clear from (2.12). For the case of 
N — 1, eq. (2.12c) is in agreement with [4]. 

Once the fermionic determinant (2.12) is positive, we see that the path-integral (2.1) is 
positive, except for the phase factor exp (iX J d 4 xX), which does not lower the ground state 
energy. This is why E(X) must have a minimum only at A = [4]. 

Before ending this section, we give the following lemma which will be of importance in 
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the next section. Note that the lagraigian £ F is rewritten as 

£f = XX - E [ +Wx Si l £&i l >) +^(C\x m J+m 3 (x'\I- J ,)+m J (e\x a J } ■ (2.13) 

Accordingly, the path-integral (2.12) is also equivlent to 

h = (n [i d x,][dx,m,mS] 



x exp 
= (detm) 2n ° 



/ d 4 xJ2{UC a Hrn + t^N)x^ + X^C^(m + t^I N ^^ 



} J det (inhN + m) 



>0 . 



(2.14) 



In other words, If in (2-12) can be computed in terms of 2-component spinors in (2.6). 
This relatinship will be helpful when we consider complicated mixed lagrangians between 
the gaugini and the quark/lepton fields in the next section. 



3. Vector-Like Theory with Supersymmetry 



We now generalize the above method to supersymmetric theories with no Yukawa cou- 
plings among quarks. Suppose we have the massive quark fermions x^ X\ C\ d 
together with the massless gluino Majorana fields A Q 7 , X^ 1 . As before, the indices 
i, j, ••• = i, 2, g = dimG are for the adjoint representations of the gauge group G. All 
the fermion-dependent terms in our lagrangian are 



+ i(T i pP. Xpi ) +i(rP« P ^ i ) +i(\* I p f, .\p I ) W(X%) +m/(C j xJ 



^ oJ X aS )-z s (X aI x & s ) 



~ KT 1 )' 



(3.1) 



The Zi and u 1 are the spin fields (squarks) in the chiral multiplets (z.,x) and 



(u\ £ l ) with 



», j, 



= i, 2, TV, in the representations N and N*. The {T 1 )^ 



are 



hermitian generators of the gauge group. As in (2.2), we can assume that is hermitian 
only with positive eigenvalues. The presence of these mixing terms with (pseudo) scalar- 
dependence have been considered to be the main obstruction for the parity-conservation in 
supersymmetric theories in the past [5] [6] [7] , because they seem to prevent us from proving 
the positive definiteness of the fermionic determinant. However, we will see that this is not 
the obstruction. After the above prescription, there is no fermion-dependent term in the 
supersymmetric lagrangian £ left over other than £ F : £ = £ F + £ B with a purely 
bosonic lagrangian £b- 
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We now consider the eigenstates for \ an d £ as in (2.4), and rewrite all the x an d 
£-depenent terms in £ F , as 



where 



titCOpp - X & \C*ew. - -^C^pJ - x ai C*W } , (3-2) 



p a . = +i{T I u*) i X a I , < = +i{z*T I ) i \J 



I — _ -IrpI \ V / 



= -i^V , uJ & = -i^A^ . (3.3) 

As usual in path-integral, we can redefine the fields in such a way that the linear terms 
in x or £ disappear. In our case, this can be done by the field redefinitions^] 

tan^tan + Waim + ifilNy 1 ]* , X aifl = X aifM + (m + ifil n)' 1 i> p a j , (3.4a) 

?<!i M = ?<Ji M + ( m + i / i/ ^) _ V^- , f.\ = x^\ + P^(m + ipI N )' l i , (3.4b) 

to have 

= - L\C^(m + iplN)/x pj , ~ ti»C«\m + ipI N )fi . jfi 

+ uJC^im + ////v) V/^. + -p. l C«\m + itiI N y\iuj.. . (3.5) 

Since m is hermitian only with positive eigenvalues, m + iplN is also diagonalizable only 
with non-zero eigenvalues, and there is no problem for defining the inverse (m + i//Jjv)~ . 
After this, £p is now 

£f = E^x,^ + E^AV + , (3.6) 

n fj. 

where £>' x £u ls the first line of (3.5) which coincides with the non-supersymmetric case 
(2.7), £\ip\ is the gluino kinetic term, while £a 2 ,^ is the A 2 and A 2 -terms after the field 
redefinition (3.4). Let us collect all of these A -dependent terms into C\. 

C x = -iA/^V -Y,K I C*M I *% J - ^A/C^M/ J A. J 

= -iX.tp^Xp 1 - \JC aP M IJ Xp J - AVC^M 7 - 7 ! . J , (3.7) 

where the matrices M M = (M/ J ), M M = (M/ J ), M = (M IJ ) and M = (M IJ ) are 
defined by 

m/ j = -( z *T^m- l T^u*) , m/ j = -(«r( / im;¥ J ) z) , 

M /J = ^M/ J , M /j = ^M/ J , m^^m + ifilN . (3.8) 



3 Note that (3.4b) is not necessarily the hermitian conjugate of (3.4a). This is related to the 
hermiticity only by the combination of + fj, and — fi. 
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Due to the antisymmetry of C af3 and C a/3 , the matrices M M , M M , M and M are all 
symmetric in i^j. Note also that 

M,,' = M j: ' = 17 „ , M t = M* = M , (3.9) 

the latter of which is confirmed by the former under J2u, which is symmetric between 
+ //<-> — /i. The M is not necessarily hermitian, and it has both real and imaginary part. 
Now C\ is rewritten as 

_ / O pj\ (M 0\ 

o H ' * s ''°(o m) ' (310b) 
o c'') = ^ c ' o c'') • (3 ' 10c) 

The original fermionic lagrangian now is Cp = £^ and the total path-integral 

to be considered in the Euclidian space is 

e- VE{x) = J [dA a In [\_dz] [dz*] [du] [du*] I F e-K^B+ax) ? (3J1) 
where the fermionic path-integral Jp is from (3.5) and (3.10) with M. replaced by M.: 
h = f [dx] W\ m m [d\] [dX\ exp [ - / d 4 x ( £ £^ + £ x ) " 

= (il/ fe] fe] fc] [4m] e^**^) / [ dA ] [rfA] e -/*** 

2 

= (detm) 2n ° JJ det (™ + ^v) [Det(^ +7W)] 1/2 . (3.12) 

At>0 

The first two factors are from the / [dx] [<^] [<£] [df ] -integral as in the non-supersymmetric 
case (2.14) now with the shifted variables x, x> an d the remaining factor is from the 

/ [dA] [dA] -integral. There is potential phase ambiguity [17] [18] for taking the square root in 
the last factor in (3.12). However, we will shortly show that there is no problem with this 
ambiguity in vector-like theories. 

Even though we can not diagonalize T/> and M. simultaneously, we still can use the 
eigenstate \v) for the eigenvalue v e 1R of the operator TJ>: 



y>\is) = +v\v) , i.e., TJ>K = 1? Q"^ = +u = i/A„ , 



(3.13) 



where the subscript v on A„ etc. denotes the eigenvalue for the four-component spinor A, 
with the adjoint index 1 omitted. As usual, we can define 

Y 5 = I 2 ® ( *° ° T ) , (3-14) 



O -Ig. 
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satisfying {r 5 , Tf> } = 0, so that an eigenstate |— v) of Tf> can be constructed by 

^[T 5 \u)]==-T 5 ^\u) = (-u)[T 5 \u)] F 5 \u)= \-u) . (3.15) 

Therefore any eigenstate for W > is always paired up with an eigenstate — v < 0. We 
now see that a conjugate state (u\ is related to \v) as follows: Consider 

A t =(^) t = (-A,, -A«) = (A,A)(_° -£ > J=A T r = JC- 1 r , (3.16a) 

^ = / 2 ®^_^ ^ , T~ X = T , |T,C]=0, [JF.M]=0 . (3.16b) 

Here A is the usual Dirac conjugate of A, and T is needed for complex-conjugation. 
Therefore 



■ f = (z/IC- 1 ^- 1 , (i/|t =TC\v) . (3.17) 



Accordingly, as in (3.15) we can confirm that 

{-v\ = (is\F 5 . (3.18) 

Other important relations needed are 

TMT- X = M* , CMC- X = M , [T 5 ,M]=0 , (3.19) 
— / — — 

as easily confirmed. Using these as well as (3.16), we get 

r 5 {ip+M)r 5 = -iTj> +M 

= + C- 1 ?- 1 ^ + M X )TC = C- X T- X { -ity + M X )TC 

= +C- 1 F- 1 (ilJ> + M) ] FC . (3.21) 

The determinant in the square root in the last factor in (3.12) can be re-expressed as the 
usual definition of the determinant in terms of exponential, trace and logarithmic functions: 

Det(i# + M) = (detM)"°(detM)™° ]J exp[(v\ In (^ + M) \v)\ . (3.22) 

As in the case of quarks, we consider only the instanton number zero background, so that 
the number of left- and right-handed zero- modes are the same: n LQ = n R(j = n , whose 
contributions in (3.22) can be computed separately, as 

(detM)"°(detMf = (det M)"°(det M)*"° = | detM | 2 ™ . (3.23) 
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As for the v ^ contributions, due to the pairing between the \u) and |— u), (3.22) is 
rewritten as 



Det {iTj> +M) 



= I det M 



i2n 



[Jexp[(i/| ln(^ +M)\v) 



u>0 



\ [ exp[(z/| In (ify + M) \v) ] 



u<0 



det M \ 2no J] exp [ (v\ In {iTj> + M) \v) ] exp [ In {%Tj> + M) \ - u) 



(3.24) 



u>0 



where the exponent in the last factor is simplified by the aid of (3.17) and (3.21) as 

{-v\ In {iTj> +M)\-v) = {u\ T 5 In {iTj> + A^)r 5 \u) 

= (H C- X T- X { In {ilj> + M)} ] FC \u) = [ (v\ In {%Tj> + M) \v) ] f 

= [{u\ ln{ip +M)\u)]* . 

Therefore (3.24) is semi-positive definite: 



Det (ip + M) = | det M \ 2n ° J] exp(i/| In (^ + A4) |i/> 



y>0 



> 



(3.25) 



(3.26) 



Combining this with (3.12), we get the semi-positive definiteness of the fermionic determi- 
nant: 

2 



det m 



|2rao 



detM 



|U0 



Q det (m + i^In) j I 

fM>0 u>0 



exp(z/| In {iTj> + At) \v) 



> . (3.27) 



Notice that the usual phase ambiguity when taking the square root [17] [18] does not 
arise here, because of the semi-positive definite expression of (3.26), as conbributions from 
| v) and | — v) always in pairs. The main ingredient in this proof is the usage of the 
eigenstate \v) with the properties of and complex conjugations, which do not require 
the diagonalization of .M, or even its commutator with Tj> . The crucial procedure we have 
relied on is the expression of the determinant in terms of exponential, trace and logarithmic 
functions, which is to be the universal definition for a determinant. 

Note that (3.27) implies only non-negativity of ip which can still be zero. This is 
because the matrix M. can depend on the scalar coordinates Z = (z, z* ,u,u*). However, 
we can further show that there exists a measurable support (a connected domain with non- 
zero measure) in the Z-space, on which If(Z) > and non-zero. In fact, consider the 
particular point Z = on which M = O due to (3.8), (3.9) and (3.10). It follows that 

2 



/ F (0) = |detm| 2no | detM 



ino 



I [ det (m + i/jlIn) \ [ 

(M>0 l/>0 



exp(z/| In (iv) \v) 



> 



(3.28) 



Once we get /f(0) > 0, then relying on the smoothness of If(Z) as a function of Z, we 
can conclude that If(Z) > on a measurable support including Z = 0. The existence 
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of a measurable support for I F > leads us to the positivity of the path-integral measure 
/ [dz][dz*][du][du*] I F > 0, and we conclude that parity is conserved in supersymmetric vector- 
like theories. 

In the above analysis, we have performed the most usual Wick rotation from the 
Minkowskian metric (+,—,—,—) into the Euclidian one (—,—,—,—) by replacing for- 
mally the coordinate x° — > ix 4 . However, this may need more care, when it comes to the 
complex conjugation of spinors. Motivated by this, we have re-confirmed our result above 
by an alternative Wick rotation into the metric (+, +, +, +). Additionally, the spinors in 
these Euclidian spaces are only 'formally' defined, in such a way that their complex conju- 
gation rule is essentially the parallel to the Minkowskian case, like the simple replacement 
x° — > ix 4 , and this is the very reason why the Feynman rules in the usual Euclidian path- 
integral is essentially the same as those in the Minkowskian. Rigorously speaking, spinors 
in the Euclidian spaces can exist only as USp(2) spinors [19], and moreover the dotted and 
undotted spinors as eigenvectors of the 75 -matrix are no longer related by complex conju- 
gations [19]. f] For this precaution, we have also reformulated the Wick rotation, such that 
the dotted X a and undotted A^ spinors in the final Euclidian space are not related to each 
other under complex conjugation, as they should be [19]. Interestingly, we have reached the 
same conclusion for the semi-positive definiteness of the determinant (3.27), even though the 
meaning of the bra- and cket- vectors are slightly modified, and all the pseudo-scalar Yukawa 
couplings with 75 acquire an extra factor of i like the AX -term in (2.1). One additional 
feature in this case we seem to rely on is that the gluino zero-modes are unstable and disap- 
pear from the physical spectrum, based on the analysis in ref. [21]. Since the details of this 
formulation is rather technical leading essentially to the same conclusion, we skip them in 
this paper. 



4. Vector-Like Theory with Broken Supersymmetry 

We mention the possibility of adding some gluino mass terms, which may be caused by 
some spontaneous, explicit, or non-perturbative breaking of supersymmetry. This is easily 
considered, by adding the gluino mass terms 



C m , X 2=m nj (\ al \ a J )+m UJ (\ ai \. J ) , (4.1) 

to our original lagrangian (3.1). Here m' = (m' IJ ) is real and symmetric. Accordingly, 
(3.10a) is now 



7 2 ®(M + m') ip®I g \ f X 



A = -(A^)pr"";;r' T ZZ" JIt ) • (4-2 
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I g h ® (M + m') J \X 



4 This situation is similar to what is called Aiyah-Ward space-time with the signature 
(+, +, -) studied in [20]. 
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This implies that the matrix M is replaced by M + m' and M by M + m'. Under this 
shift, the relations such as = M are intact. Eventually (3.27) is now replaced by 



4 



det m 



1 2no I 



detM 



ino 



Y[ det (m + ifiI N ) JJ exp(z/| In + X') 

/K>0 i^>0 



> , (4.3) 



where M.' is a 4g x 4g matrix similar to M. defined by 

'm' O \ 
O m'J ' 



A4' = A4 + J 2 



M' ] =M' 



(4.4) 



Hence the presence of m' does not affect the semi-positive definiteness of the fermionic 
determinant. Accordingly, we can also show that J F > on a measurable support in the 
Z -space, and therefore we conclude that parity is conserved also in broken supersymmetric 
vector-like theories with non-zero gluino masses. 



5. Concluding Remarks 

In this paper we have shown the conservation of parity in supersymmetric vector-like the- 
ories. The main body of our proof is the confirmation that the determinant as the fermionic 
path-integral Ip is positive and non-zero on a measurable support under the bosonic inte- 
gral / [dz][dz*~\[du][du*~\. We have also seen that the quark-gluino-squark mixing terms in the 
supersymmetric theory pose no problem. The supersymmetric vector-like theory seems to 
avoid the problem with these mixing terms thanks to parings between the eigenstates of the 
Dirac operator, despite of the complication caused by the mixing with gluini. 

In our analysis of the fermionic path-integral, we first integrated over the quark fields 
X) Xi £) making the computation more organized, instead of integrating over the gluino 
field first as in [6]. By so doing, we have seen that the final gluino path-integral is less involved 
and more controllable, in particular when we need to consider the Majorana gluino determi- 
nant which used to have subtlety with 75 -pseudo-scalar couplings. We have understood that 
the usual phase ambiguity in the square root of the determinants for a 2-component spinor 
[17] [18] does not arise in a vector-like theory, due to the pairing between the integrals over 
dotted and undotted spinors, combined with the pairing between the positive and negative 
eigenstates of the Dirac operator, including the zero-modes. 

We have also studied the effect of the gluino masses, as a result of either spontaneous, 
explicit, or non-perturbative breaking of supersymmetry, and reached the conclusion that 
parity is also conserved in these cases with broken supersymmetry. In principle, we can 
also consider the squark masses caused by the supersymmetry breaking, but these terms are 
purely bosonic affecting only Cb, so that they are not expected to change our analysis or 
result in this paper. 
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Note that our result relies on the original method in [4], namely we deal only with bosonic 
parity-breaking observables in path-integrals. Therefore our result does not cover the parity- 
breaking via fermionic observables [12] which do not acquire the imaginary unit % under 
the Wick-rotation. 

Our result here seems to contradict with ref. [9], which shows that the vacuum structure 
is disturbed, when there are more colors than flavors: Nf < N. However, there is actually 
no conflict, because we interpret this as the result of the masslessness of quarks treated 
in ref. [9]. In our system, due to the massive quarks from the outset, the Witten index 
Tr (— 1) F is well-defined and non-zero [10], e.g., Tr (— 1) F = N for the SU(N) gauge 
group. Hence we expect no breaking of supersymmetry even at the non-perturbative level, 
which would have ruined the foundation of our proof. Since the topological stability due 
to the well-defined Tr( — 1) F is reliable for massive quarks, it is quite natural that chiral 
symmetry or parity symmetry is also conserved. Additionally, our parity-conservation is also 
consistent with the result of [11] about the phase with unbroken discrete axial symmetry. 
Some subtlety arises, when the masses of the quarks become zero, because in such a case 
the Witten index Tr(— 1) F is no longer well-defined [10] [22], and therefore the vacuum 
loses its stability against chiral or parity breakings as in [9]. From this viewpoint, we see no 
contradiction of our result with refs. [9] [22], in which the masses of quarks are zero from the 
outset. 

We have seen that the supersymmetry breaking, if its only effect on fermions is the gluino 
masses, does not alter the conservation of parity. Even though this statement seems contra- 
dictory with the previous paragraph, we understand that the breaking of supersymmetry lifts 
the vacuum energy higher than the original supersymmetric and parity-conserving vacuum 
with E{0) = 0. 

We can try to apply our method to other arbitrary global or discrete symmetries, such as 
baryon number, in a supersymmetric vector-like theory, using the prescription using upper 
bounds for fermionic propagators in ref. [23]. However, there seems to be an obstruction 
caused by the zero-ness of mass eigenvalues, e.g., our matrix M. in (3.10) hitting zeros, 
that upsets the upper bound for fermionic propagators, undermining the foundation for the 
stability of fermion-anti-fermions Greens functions against symmetry-breaking parameters 
[5] [23]. In other words, our method using the semi-positive definite fermionic determinant 
is powerful only for parity symmetry, or other symmetries based only on the vacuum to 
vacuum amplitude. 

We are grateful for J.C Pati for bringing about the problem with parity conservation in 
supersymmetric QCD, and for other helpful discussions. We are also indebted to M. Luty, 
R. Mohapatra, Y. Shamir, and E. Witten for important suggestions and comments. Addi- 
tional acknowledgement is due to M. Cvetic for a copy of reference [6], and other important 
discussions. Special acknowledgement is for C. Vafa for pointing out mistakes in an earlier 
version of the manuscript. 
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